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<N : ABSTRACT 



We report the first detection of [C I] ( 3 Pi — > 3 Pq) emission at 609 /xm 



> 

in 

in a region of the Small Magellanic Cloud (N27). Environments poor in 
heavy elements and dust such as the SMC are thought to be dominated 

O 

O - 

submillimeter neutral carbon emission has been detected. Studying the I[ci]/Ico 
ratio in several sources spanning more than an order of magnitude in metallicity, 
Z, we find that the Ifcrj/Ico ratio increases for decreasing Z. The existence of 
such a trend points to a photodissociation origin for most of neutral carbon in 
molecular clouds, in agreement with standard PDR models. We also report 
ISO FIR spectroscopic observations of N27, and use them to derive its physical 



properties. Comparison between the density and radiation field revealed by 
FIR diagnostics (n ~ 300 — 1000 cm -3 , % uv ~ 30 — 100xo) an d those derived 



from millimeter and submillimeter data (n ~ 10 5 cm -3 , % uv ^ 30xo) suggests 
that the FIR lines originate in more diffuse gas, and are perhaps dominated 
by the interclump medium. Regardless of the cause, analysis of the FIR and 
mm-submm data produces a discrepancy of two orders of magnitude for the 
density of this source. 
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1. Introduction 

Studies of the interstellar medium (ISM) in the Small Magellanic Cloud (SMC) allow 
us to investigate star forming clouds that are poor in heavy elements. The SMC is indeed 
an excellent choice to study this phenomenon, because it is nearby (D ~ 61 kpc; Laney 
Sz Stobie 1994), has an unobscured line of sight and very small internal obscuration, and 
possesses the lowest metallicity among the nearby members of the Local Group (Z ~ Z /1O; 
Dufour 1984). Because of their low metallicity, morphological irregularity, and active star 
formation, the Magellanic Clouds resemble primeval galaxies. Therefore, studies such as 
this will eventually allow us to draw conclusions about those much more elusive targets: 
the high-z primitive galaxies. 

The depletion of heavy elements, and the consequent low dust-to-gas ratio 
(.D/Gsmc ~ 1/8-D/Gmw; Bouchet et al. 1985), affects the ISM by enhancing 
photodissociation. Because of its high abundance (i.e., strong self-shielding) and cross- 
shielding with atomic hydrogen lines, molecular hydrogen is mostly unaffected by Z. The 
usual tracer molecules such as CO, however, are efficiently destroyed in dust-poor systems. 
A much larger fraction of the ISM in these systems is thus expected to be occupied by 
photodissociation regions (PDRs; e.g., Tielens & Hollenbach 1985; Sternberg & Dalgarno 
1995). Consequently, we expect the intensities of the spectral lines associated with the 
PDR (e.g., [C II], [C I]) to be enhanced when compared to molecular lines (e.g., 12 CO 
( J — 1 — > 0)), resulting in increasing I[cii]/Ico and I[ci]/Ico ratios (Bolatto, Jackson, & 
Ingalls 1999). 

The immediate purpose of this work is twofold: 1) to test for the existence of an 
increasing I[ci]/Ico ratio with decreasing metallicity Z (§ |3.1| ), and 2) to investigate the 
power of far-infrared (FIR) spectroscopic diagnostics and compare their results to those 
obtained using millimeter and submillimeter lines (§|3.2|). An I[ci]/Ico dependence on Z 
could be useful for determining metallicities in high-z sources using only radio observations. 
The trend in Ifcrj/Ico with Z is also a test for the origin of neutral carbon inside molecular 
clouds. An increasing I[ci]/Ico ratio for decreasing Z points to a PDR origin, while a flat 
I[ci]/Ico ratio suggests that C° is produced inside the CO cores by means other than UV 
photons. 

Our SMC source, the N27 nebula, was first cataloged by Henize (1956). A later survey 
by Davies, Elliott and Meaburn (1976) calls it DEM40. Its associated IR peak is also 
known as LIRS49 (Schwering & Israel 1989). This region was mapped by the Swedish-ESO 
Submillimeter Telescope (SEST) in the ( J — 1 0) and ( J — 2 — > 1) transitions of CO 
(Rubio et al. 1993), and studied in detail in several molecular transitions by Heikkila et 
al. (1999). Albeit a weak emitter when compared to Galactic sources, the molecular cloud 
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associated with N27 features some of the strongest CO ( J = 1 — > 0) emission in the SMC 
(Israel et al. 1993). 

2. Observations 

2.1. Far Infrared Data 

We observed N27 using the Long Wavelength Spectrometer (LWS) onboard the Infrared 
Space Observatory satellite (ISO), at R = A/AA ~ 200. These LWS data are a 4160 
second-long full-grating scan obtained on 20 December 1997 as part of the LOMETPDR 
project. The instrumental characteristics are described by Clegg et al. (1996), and the 
calibration of the ISO LWS instrument is described by Swinyard et al. (1996). The data 
pipeline was version 8.7 and the observations were processed using ISAP version 2.0. We 
averaged the spectra across scans and scan direction, but not across detectors, with the 
default wavelength spacing of 0.135 /im. The averaging was done after 2.5a clipping of 
outlying points. For the purposes of obtaining an I FIR estimate only, we joined the detectors 
together by offsetting them with the dark correction method in ISAP, using detector LW1 
as the reference. Table [I] summarizes the spectroscopic results. 

We can assess the calibration of the LWS data by comparing it with previous [C II] 
measurements of N27 featuring a similar beam size (Israel & Maloney 1993). The intensity 
measured by FIFI at the LWS pointing is I[cn]~ 9 x 10~ 4 erg s" 1 cm" 2 sr" 1 . The LWS 
measured intensity after applying the recommended correction factors for beam size and 
extended source is I[cii]~ 9.5 x 10 -4 erg s _1 cm -2 sr _1 , in remarkable agreement with the 
FIFI observations. The measurements listed in Table [I] include both corrections. 

2.2. Submillimeter Data 

We observed the [C I] ( 3 Pi ^ 3 Po) transition of neutral atomic carbon at v ~ 492.1607 
GHz (609 /jiii) using AST/RO, the Antarctic Submillimeter Telescope and Remote 
Observatory located at the Amundsen-Scott South Pole base (Stark et al. 1997). The 
observations were obtained on 3 to 7 July of 1998, using the AST/RO SIS quasi-optical 
receiver, with a system temperature T sys ~ 1800 — 2400 K. The backend was the 2048 
channel low resolution AOS. The spectra were observed in position switching mode, 
chopping 15' in Azimuth (same as R.A. at the pole). At 492 GHz the telescope beam had 
HPBW~ 3.8' ± 0.3. The forward efficiency determined from skydips was 70%, and it is 
assumed to be identical to r] mb . The pointing was better than 30". A 6' x 6' region centered 
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on this source was mapped on a 30" grid, considerably oversampling the beam. The data 
were calibrated using the standard procedure for AST/RO, which includes sky, hot and cold 
load measurements every 30 minutes, and processed using the COMB package. The total 
accumulated integration time, after rejecting bad spectra, was 63 hours. The spectrum 
obtained after averaging all the observations together is shown in Fig. [l]. A gaussian fit to 
the data gives v cen = 115.2 ± 1.4 km s _1 , Av = 9.3 ± 3.3 km s _1 , and T max = 14 ± 4 mK 
(main beam brightness). The location of the [C I] line agrees very well with the strongest 
CO component, for which v cen = 114.5 km s _1 (Rubio et al. 1993). 

The relatively weak [C I] emission from N27, combined with its small size and 
consequent beam dilution, produced a low signal-to-noise map. Thus the analysis in § |3~T 
will be focused on the averaged spectrum for the region. 



3. Results and Discussion 

3.1. Metallicity and the I[cii/Ico Ratio 

The prime motivation for this study is to determine the I[ci]/Ico intensity ratio in 
the SMC, an environment that because of its low metallicity and low dust-to-gas ratio 
should be dominated by PDRs (cf., §|l|). To do this we compare our [C I] data with the 
12 CO (J = 1 — > 0) map by Rubio et al. (1993). This map covers only a relatively small 
fraction of the area mapped in [C I], but it probably contains most of the CO emission: to 
obtain a I[ci]/Ico ratio we assume that there is no CO emission outside this region. The 
CO emission associated with N27 has two components, at 114.5 and 126.9 km s _1 . For this 
analysis we have integrated the CO and [C I] data over the velocity range corresponding 
to the strongest component, at ~ 114 km s _1 . In order to compute a ratio we have 
convolved the CO ( J = 1 — ► 0) integrated intensity data to the angular resolution of the 
[C I] map. Subsequently we sampled the convolved CO map at the positions observed in 
[C I] , excluding those positions where [C I] spectra were discarded. These samples were 
averaged with the same weighting scheme used for [C I]. This procedure allows us to obtain 
a reliable average ratio for the N27 region, relatively insensitive to uncertainties in the 
[C I] beam size and pointing. The resulting intensities, integrated over the velocity interval 
107—119 km s _1 , are I [CI] ~ 0.132 ± 0.027 and I co ~ 0.517 ± 0.052 K km s _1 , yielding 
a ratio I[ci]/Ico ~ 0.26 ± 0.06. The corresponding cooling ratio (the ratio of intensities 
in erg s" 1 cm -2 sr _1 ) is Ifcii/Ico ~ 19.9 ± 4.5. The uncertainties are la and include the 
statistical errors as well as a 10% la calibration uncertainty in the intensities, added in 
quadrature. 
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Figure |2] compares the I[ci]/Ico ratio in N27 to a sample of molecular clouds associated 
with star-forming regions in systems of different metallicities. This plot suggests that 
there is an increasing trend in the I[ci]/Ico ratio with decreasing metallicity Z, such that 
I[ci]/Ico oc Z~ ' 5 . The scatter around this line is relatively large, as shown by the data for 
the Large Magellanic Cloud. Theoretical expectations notwithstanding, this scatter has 
made it difficult to identify a trend in past studies over a narrower range of metallicities 
(e.g., Bolatto et al. 2000a). 

Shown also in Fig. |2] is a model for carbon in unresolved PDRs (Bolatto et al. 1999). 
This model consists of intensity ratios computed over a distribution of clumps with uniform 
excitation, where each spherical clump is divided in three concentric regions of C + , C°, 
and CO. The sizes of these regions depend on the metallicity (or the dust-to-gas ratio) of 
the ISM. While necessarily crude this model incorporates the effects of clumping, and the 
increase in the physical size of the PDR for decreasing metallicity. The growth of the PDR 
is mostly caused by the diminishing attenuation of the UV radiation, as the dust-to-gas 
ratio of the ISM decreases for smaller metallicities. In this scenario, the primary reason 
why the Irci]/Ico ratio increases for decreasing metallicity is the reduction and eventual 
disappearance of the CO cores in the clumps. 

Are the models for the Ircii/Ico ratio as a function of Z quantitatively consistent with 
the measurements? Bolatto et al. explore two possibilities: C° produced in the growing 
PDR (model A), and C° restricted to a region whose size does not change with metallicity 
(model B). This latter model predicts an essentially flat I[ci]/Ico ratio (I[ci]/Ico oc Z +0A ). 
A similar result (constant I[ci]/Ico) would be obtained if C° was produced inside the 
CO cores via unidentified chemical reactions. The slope predicted by model A in this 
metallicity range is slightly steeper than that suggested by the data, I[ci]/Ico oc Z~ - 8 . 
Whether this difference is actually significant is difficult to assert given the paucity of 
observations and the intrinsic scatter in the data. Nevertheless, several possibilities could 
explain this discrepancy: 1) an increase in CO excitation for decreasing metallicity (i.e., 
hotter molecular gas for diminishing Z), 2) clump-to-clump shielding preventing the 
efficient photodissociation of CO, 3) a dust-to-gas ratio not proportional to the oxygen 
abundance (but this would require the dust abundance to be a slower function of Z, e.g. 
D/G~ Z ' 5 , when observations suggest D/G~ Z 2 \ Lisenfeld & Ferrara 1998), and 4) some 
fraction (~ 30%) of the C° produced inside the CO cores instead of the PDR. We conclude 
that, within the bounds imposed by the limitations of the models and the paucity of the 
measurements, the observed trend for the I[ci]/Ico ratio with metallicity shown in Fig. |2| 
agrees with theoretical expectations for mostly photodissociation-produced C°. 
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3.2. Physical Conditions 

The FIR spectroscopic data can be used to determine the conditions prevalent in N27. 
Using Table [TJ and the PDR models by Kaufman et al. (1999) we will find the UV field, 
Xuv, and the gas density, n. The important FIR ISM diagnostics computed by Kaufman 
et al. are: 1) the [O I] ( 3 P ^ 3 Pi)/( 3 Pi ^ 3 P 2 ) line ratio (Wfe ~ 0.051 ± 0.006, 
according to our measurements), 2) the [O I] ( 3 Pi — > 3 P2) to [C II] ( 2 P3/2 - ^Pi/2) line ratio 

(Wliss ~ °- 42 ± °- 01 )' 3 ) the ratio of [° A ( 3p i ^ 3p 2) + [C II] ( 2 P 3/2 ^ 2 Pi/ 2 ) to the FIR 
continuum (using the LWS spectrum we estimate Ifir ~ 1.1 X 10~ 2 erg s" 1 cm" 2 sr" 1 , 
yielding (I 63 + Ii 58 )/I F ir ~ 0.012 ± 0.002), and 4) the [C II] ( 2 P 

3/2 — >2 Pi/2) intensity 

itself. These joint four diagnostics agree very well, placing N27 in the low density, 
n ~ 300 — 1000 cm -3 , moderate UV field, Xuv ~ 30 — 100xo, region of the n — % uv diagrams 
(Xo is the interstellar UV field in the vicinity of the Sun, xo — 1-6 X 10~ 3 erg cm -2 s _1 ; 
Habing 1968). The corresponding PDR temperature in these conditions is T p d r ~ 200 K. 
The FIR continuum dust temperature, derived from the LWS spectrum, is Td ~ 40 K. 

The density and radiation field derived using millimeter and submillimeter lines, 
however, strongly disagree with the FIR results. This lack of agreement is not surprising, 
perhaps, since these lines probably originate in different regions of the ISM. A millimeter 
wavelength multiline excitation analysis using several molecular species was performed by 
Heikkila et al. (1999), yielding n ~ 10 5 cm -3 , T ^ 15 K, and visual extinction Ay ~ 1 for 
N27. The I[ci]/Ico ratio determined in §[3j] strongly constrains the density according to the 
calculations by Kaufman et al. (1999), also resulting in n ~ 10 4 — 10 5 cm -3 (albeit over a 
much larger area). Finally, using the 12 CO ( J = 1 — > 0) observations by Rubio et al. (1993) 
convolved to the FIFI beam we compute an accurate I[cii]/Ico ratio at the LWS pointing, 
yielding I[cii]/Ico ~ 7000. This ratio alone does not constrain either parameter very well, 
but implies Xuv > 100xo arid when used in combination with the measured I[ci]/Ico ratio 
results in % uv ~ 10 4 — 3 x 10 5 Xo- This is not a fair comparison, however, since both the 
I[cii]/Ico and I[ci]/Ico ratios will be affected by the low metallicity and extinction of this 
source. Fortunately, Kaufman et al. (1999) provide some calculations for small Z and Ay- 
Using their model tracks for Z = 0.1 and Ay = 1 we find n ~ 10 5 cm" 3 and Xuv ~ 30xo 
for N27, based on I[cii]/Ifir ~ 10~ 2 and Ico/Ifir ~ 10~ 6 . These measurements are also 
consistent with nco > n c+ at a lower radiation field, if we allow the model tracks to be 
shifted in the abscissa as described by the model authors. 

At any rate, the diagnostics including [C I] or molecular line observations are 
inconsistent with the density n ~ 10 3 cm" 3 implied by the FIR transitions. Because in 
general the FIR lines are generated closer to the surface of the cloud, this discrepancy may 
indicate a real density gradient in the gas. Another possibility is that the mm-submm and 
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FIR lines originate in different regions of the cloud, with the latter dominated by material 
in a more diffuse interclump medium, or perhaps associated with the H II region rather than 
the molecular cloud. A third possibility is that the FIR line ratios are significantly affected 
by Z or Ay, and using model results computed for Galactic metallicity and large Ay is 
inappropriate. If this is the cause of the density discrepancy, however, it is remarkable that 
all four FIR diagnostics agree on essentially the same value for the density and radiation 
field. Independent of which explanation is correct, we conclude that: 1) the FIR and 
mm-submm measurements roughly agree on the radiation field intensity, although it is 
extremely important to use models with the appropriate Z and Ay, and 2) the mm-submm 
data systematically indicate a density about two orders of magnitude higher than that 
derived from FIR fine structure data. 



4. Summary and Conclusions 

We report new FIR and submillimeter observations of the star forming region N27 in 
the Small Magellanic Cloud. Using our measured intensities for the FIR continuum and 
[O I] and [C II] fine structure transitions, and comparing with existing model calculations 
(Kaufman et al. 1999), we derive the density and radiation field in N27 to be n ~ 300 — 
1000 cm~ 3 and Xuv ~ 30 — lOO^o- in contrast, the density derived using millimeter and 
submillimeter data is two orders of magnitude higher, n ~ 10 5 cm -3 , while the radiation 
field is similar or lower, % uv ^ 30xo- The cause for this discrepancy is not clear, but we 
suggest that the FIR and the mm-submm lines arise in different parcels of gas, with the 
former perhaps dominated by the interclump medium or the H II region. 

We have detected [C I] ( 3 Pi — > 3 Po) emission from the Small Magellanic Cloud for 
the first time. This is the lowest metallicity source in which submillimeter neutral carbon 
emission has been detected. Low metallicity, dust-poor environments are thought to be 
dominated by PDRs. Despite the factor of ~ 5 in metallicity spanned by N27 and Orion, 
however, the I[ci]/Ico ratio in N27 is only modestly larger than the I[ci]/Ico ratio in 
Orion. Nevertheless, when we study the ratio in several sources spanning more than an 
order of magnitude in metallicity, we find a noisy but convincing increasing trend in the 
I[ci]/Ico ratio for decreasing metallicity (Fig. g). This trend is somewhat flatter than, but 
in rough agreement with, the modelling results by Bolatto et al. (1999). The existence of 
an increasing trend points to a photodissociation origin for most of neutral carbon inside 
molecular clouds, in agreement with standard PDR models. 
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Table 1. ISO LWS Line Intensities in N27 



Specie 


Transition 


Wavelength 


Beam Size 


ESCF a 


Line Flux 


Continuum 


Line/Cont. b 








(jaa) 


(sr) 




(erg s^ 1 cm -2 sr -1 ) 


(MJy sr" 1 ) 




[O III] 


( 3 P 2 - 


> 3 Pl) 


51.815 


ii.8 x icr 8 


0.85 


1.34±0.22 x 10~ 5 


191 


0.20±0.03 


[O I] 


( 3 Pl- 


> 3 P 2 ) 


63.184 


13.0 x icr 8 


0.83 


3.93±0.09 x 10- 5 


237 


0.72±0.02 


[0 III] 


( 3 Pl- 


> 3 Po) 


88.356 


11.8 x 10~ 8 


0.67 


9.64±0.85 x 10~ 6 


293 


0.13±0.01 


[O I] 


( 3 Po- 


» 3 Pl) 


145.525 


6.7 x icr 8 


0.58 


2.01±0.23 x 10~ 6 


415 


0.06±0.01 


[C II] 


( 2 P 3 /2 


- 2 Pl/2) 


157.741 


6.7 x icr 8 


0.58 


9.46±0.08 x 10~ 5 


394 


3.14±0.04 



a Rccommcndcd extended source correction factor (ISO LWS Handbook online vl.0). 

b Line to continuum ratio at line center (unresolved lines, AA = 0.29 (im for A < 80 (an. Otherwise AA = 0.60 ^m). 
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Fig. 1. — Averaged [C I] spectrum of N27 obtained by AST/RO (Hanning-smoothed to 
Av ~ 2.1 km s -1 ). Superimposed are the rescaled Gaussian 12 CO components measured by 
Israel et al. (1993). There is a hint of the 127 km s _1 component in the [C I]. 
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Fig. 2. — The [C l]/CO line intensity ratio as a function of metallicity, for molecular clouds in 
the Galactic Center (GC), the Orion region, the LMC, IC 10, and the SMC. The continuous 
line is a linear fit to the data: log[I [C i]/Ico] = -0.47 * (12 + log[0/H]) + 5.0. The dashed 
line is Bolatto et al. (1999) model A (PDR-produced neutral carbon): log[I[ C i]/Ico] = 
—0.82 * (12 + log[0/H]) + 8.2. The references for metallicity and I[ci]/Ico ratio for each 
source indicated in the figure are: (1) Afnerbach et al. 1997, (2) Ojha et al. 2000, (3) 
Tauber et al. 1995, (4) Dufour 1984, (5) Bolatto et al. 2000b, (6) Lequeux et al. 1979, (7) 
Bolatto et al. 2000a, and (8) this work. 



